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Abstract 

We examine cross-correlations (CCs) in the tunneling currents through two parallel interacting 
quantum dots coupled to four independent ferromagnetic electrodes. We find that when either one 
of the two circuits is in the parallel configuration with sufficiently strong polarization strength, a 
new mechanism of dynamical spin blockade, i.e., a spin-dependent bunching of tunneling events, 
governs transport through the system together with the inter-dot Coulomb interaction, leading to 
a sign-reversal of the zero-frequency current CC in the dynamical channel blockade regime, and to 
enhancement of positive current CC in the dynamical channel anti-blockade regimes, in contrast 
to the corresponding results for the case of paramagnetic leads. 

PACS numbers: 72.70.+m, 73.23.Hk, 73.63.-b, 72.25.Rb 
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Quantum noise cross-correlation (CC) in multiterminal mesoscopic devices far from equi- 
librium has recently become an active issue because it characterizes the degree of correlation 
and the statistics of the charge carriers.- 1 ^ It is believed that in a noninteracting system, the 
current CC between different normal-metallic leads is always negative due to the fermionic 
statistics of electrons; 3 this has been confirmed experimentally in a Hanbury Brown-Twiss 
setup.- On the other hand, much theoretical work has predicted that current CC may become 
positive in the following situations: for a hybrid superconductor-normal system; 5 for a sys- 
tem with three-terminal ferromagnetic leads;- also for a system in the Coulomb-interaction 
mediated regime;- and as a result of feedback effects of external voltage fluctuations. 8 

Very recently, Coulomb-interaction-induced positive current CC was experimentally ob- 
served in a double quantum dot (QD) system, with two parallel QDs coupled via an inter-dot 
Coulomb interaction, and also coupled to four independent electrodes and two independent 
gates, as shown in Fig. I. 9 It was found that when the system is biased on the specific 
occupation parameters, (1,0) and (0,1) [(M,N) denotes the electron numbers in the top 
and bottom QDs], by means of tuning two gate voltages, the electrons occupying one dot 
during sequential tunneling through that QD can enhance either the tunneling-in rate or 
the tunneling-out rate of the other dot due to the Coulomb interaction (dynamical Coulomb 
anti-blockade), resulting in the appearance of a positive CC between the currents through 
the two QDs (i.e. electron bunching).-^ However, the current CC still remains negative 
in the transport regimes, (0,0) and (1,1), due to the dynamical Coulomb blockade effect 
(electron anti-bunching) . 

In this Letter, we analyze the current CC for the same structure as the one studied in 
Refs. , but with ferromagnetic leads. Experiments with a strongly interacting QD and 
ferromagnetic contacts have recently been performed successfully,— and it is not difficult to 
make such a structure with parallel-coupled QDs connected to four independent ferromag- 
netic leads, as shown in Fig. 1. Our main finding is that sufficiently polarized contacts can 
lead to positive current CC even in the Coulomb blockaded region (0,0) if only one of the 
two circuits is in the parallel configuration. Similar to results for the three-terminal QD 
with ferromagnetic electrodes studied by Cottet et a/,- the present result also stems from 
dynamical spin blockade, associated with dynamical channel blockade. Furthermore, we find 
that the sign of current CC reduces to the result for paramagnetic electrodes with increasing 
spin relaxation rate. 
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Figure 1: (Color online) Schematic diagram of the system. 



The system we study here consists of two parallel single-level QDs [top and bottom 
(i = t, b)} with energies and inter-dot Coulomb interaction U, which are coupled to four 
collinearly spin-polarized leads with net spin-independent tunneling rates 7$ and polarization 
strengths Pi (0 < Pi < 1) (Fig. 1). Accordingly, the ferromagnetism of the leads in the 
top/bottom circuit can be accounted for by the spin- dependent tunneling rates: Fl^ = 
r^T = 7»(1 + Pi) and F Lil = T m = 7^(1 - p^ for the parallel (P) configuration; T m = 
T Ri i = 7j(l + p^ and T Li ^ = T Ri ^ = 7^(1 — Pi) for the anti-parallel (AP) configuration. 
Moreover, a constant spin-flip scattering rate, j s f, is introduced to model spin relaxation. 
We also assume infinite on-site Coulomb repulsion to guarantee no double occupation on 
each dot. Therefore, there is a total of 9 states in this system: no electron in the two QDs, 
1 00), the top QD occupied by one electron with spin a, \a0), the bottom QD occupied by 
one electron with spin a, [Oct) , and either QD occupied by one electron with spin a (a'), 
In the sequential-tunneling limit, electronic transport through such double QDs can 
be described by the rate equation:-*^ 



in which p a p is the occupation probability for the state \a(3). The term M a ^^ = 
^ Q /0/ M a ip^ a p gives the total loss rate for the state \a/3), while the term M a ^^pi (a[3 7^ a' (3') 
gives the total rate for transitions between the two states \a/3) and \a'(3'). All these 
terms depend on the spin-flip transition rate and the spin-dependent rates for tunneling- 
in and tunneling-out between a QD % and its lead rji (77 = L, R) when the other QD is 
either empty or occupied, rj ff = T vi(T f^) and fJ CT = T^f^Si + U), where /+(e) = 




(1) 
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{1 + exp[(e — [irp) I ksT] } 1 is the Fermi function of lead rji with chemical potential fi V i and 
temperature T, and /~(ej) = 1 — /i(ej). For example, the total loss rates for the states \a0) 
and \aa') are M ct0 , ct0 = - E^Q^a+f J T +f Jj-T-/ and M CT(T>CT , = - £„(f ~ 6ct ,)- 2 7 s /, 
respectively. Here, a symmetric bias voltage is assumed to be applied between the left 
and right leads: Hu = —Pm = el^/2. 

The steady-state value of p°^ under finite bias voltage Vi is obtained by solving dp a p/dt = 
in Eq. ([I]). Then the spin-dependent currents flowing through the right leads can be 
evaluated as (we use H = e = ks = 1) 

iti(t) = r^ T p T0 - r^ T poo + r raT (p T T + Pn) 

-frar(PoT + Poi), (2) 

-fWPto + P|o), (3) 
and and are obtained by interch anging in Eqs. ([2]) and ([3]), respectively Using 



the techniques developed in Refs. 



13 



14j, we introduce spin-resolved current matrices 



for the two circuits according to Eqs. ([2]) and ([3]), and then apply them to the steady-state 
solutions of Eq. (CD) to calculate the average currents I t (b) = E CT h{b)a an d the spin- dependent 
noise CC S ta ,ba'(uj) defined as 

/+oo 
dre^(A/ tCT (r)AJ 6a ,(0)), (4) 
-oo 

with Al ia (t) = Ii a (t) — (Iia). The total zero- frequency current CC is S t b = Eo-o-' ^to-M' (®) ■ 

In the following calculation, we focus on the situation of two identical circuits, Pt = Pb = P 
and 7* = 76 = 7. The other parameters are: U = 67, T = 7/2 and a given bias voltage 
V t = Vb = V = IO7 (hereafter, in the text and in the figures we use 7 as the energy unit). 
Figure 2 exhibits the calculated zero-frequency current CC dependence on the energy levels 
of the two QDs, e t and for paramagnetic leads (a), and for ferromagnetic leads with 
strong polarization, p = 0.6, and no spin-flip scattering 7 s j = (b-c). In the paramagnetic 
case, the current CC is, as expected, negative for the occupation regimes, (0,0) and (1,1), 
while it is positive for the regimes, (1,0) and (0,1). It is already known that this positive 
CC stems from interaction-induced enhancement of the tunneling rate.- To reveal this with 
greater clarity, we have derived an analytic expression for the current CC for the regime 
(1,0) (for example, setting e t = —5.5 and Sb = 0) by assuming an enhanced tunneling-in 
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Figure 2: (Color online) Zero-frequency cross-correlation dependence on both gate voltages for 
transport bias voltage V = IO7. (a) is plotted for paramagnetic leads; (b-c) are for ferromagnetic 
leads with p = 0.6 and without spin-flip scattering 7 s j = 0. Both circuits are in the P configuration 
(b), the AP configuration (c), and one circuit is in the P configuration and the other one is in the 
AP configuration (d). 



rate for the bottom QD under the condition that the top QD is occupied by an electron, 



Lba 



x (0 < x < 1/2) [correspondingly, the tunneling-out rate is T Lba — 1 — x; the other 
rates are either 1 or 0]. In the limit of zero temperature, the analytic result is given by 

2 89 x 2 + 2813 x- 1890 n 

S '^-Ti-< X (^Tsf >0 ' (5) 

Surprisingly, we find that the sign of current CC for the transport regime (0,0) is changed 
in the presence of ferromagnetic electrodes [Fig. 2(b-c)], and its largest values are obtained 
in the case when both two circuits are in the P configuration (P-P configuration). Fur- 
thermore, the positive CC noise for the transport regimes, (1,0) and (0,1), is enhanced in 
comparison with those in the case of paramagnetic leads, if either of the two circuits is in 
the P configuration. This sign-reversal and the enhancement can be explained in terms of 
the associated effect of dynamical channel blockade and a mechanism of dynamical spin 
blockade, as suggested by Cottet et al.. 6 



Let us focus on the P-P configuration. Here, we assume that up spins are in the majority. 
Correspondingly, the dwell time of down spins on the dot is longer than that of up spins 
since the down spin tunneling rates (slow electrons) are much lower than those of up spins 
(fast electrons). Considering further the fact that the total number of electrons with up spin, 
rijf , is equal to that of electrons with down spin, n^, in this configuration, we can infer that 
every tunneling event of a down-spin electron flowing through a QD must follow several 
consecutive tunneling events of up-spin electrons flowing through the QD, i.e. up spins 
bunching on both circuits. This indicates that for most time intervals, the sign of AI t ^(r) is 
the same as that of Ai^(O), but is opposite to that of Aibj(O), which is responsible for the 
occurrence of positive CC between up spins, S^^O); but negative CC involves down spins, 
SVf,6.i.(0), for the transport regimes, (0,0), (1,0), and (0,1). We have obtained analytical 
expressions for the various spin-resolved CCs as: 

WO) = ~— 7(l+p) 2 7(1 ~ 5P | + 27 2 t (6) 
Q fn\ 2 M ^ 2 7(l + 5p) + 2 7s/ 

StlMK ) = — 7777(1 — p) ■ -i-, (7) 



125 ,v ^ 2 7s/ + 7 (l- F 

WO) = Wo) = ^(i -p ) 27j/ + 7(1 _ p2) - 



(8) 



and for the total CC as: 



8 _ 7 (1 - hp 2 ) + 2 lsf 



in the transport regime (0,0) at zero temperature. Clearly, the positive sign of St&(0) stems 
from the up- up CC S^^O) > for sufficiently strong polarization p > l/y/5 ~ 0.45 
in absence of spin-flip scattering 7s / = (note that the other spin-resolved CCs are all 
negative). Moreover, retaining spin coherence during tunneling is a necessary condition for 
such bunching of up spins. We exhibit the effect of spin-flip scattering on the zero-frequency 
current CCs in Fig. 3. It should be noted that spin-flip scattering strongly influences the 
up-up CC once 7 s j is comparable to the tunneling rate 7; and in the limit of infinitely high 
7 S /, the CCs tend to the corresponding values of the paramagnetic case for any value of 
the polarization and any configuration [see Eq. (jHJ)]. In particular, 5#,(0) becomes negative 
around 7s / ~ 2 7 in Fig. 3(a). Moreover, a weak polarization, p, can not cause sufficiently 
strong up-up spin bunching to overcome the negative CCs, S^^O) and ^^^(O), leading to 
a small negative value of ^(O), as shown in Fig. 3(d) for the results with p = 0.3. 
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Figure 3: Spin-resolved noise cross-correlations as a function of spin-flip scattering, j s f, for the case 
of ferromagnetic leads in the P-P configuration, (a-c) are plotted for p = 0.6 (solid lines), p = 
(thin line), and the chosen values of energy levels of the two QDs corresponding to the various 
transport regimes, (0,0), (1,0), and (1,1), respectively; (d) is plotted for p = 0.3 in the transport 
regime (0,0). The Inset in (a) shows the corresponding results for the transport regime (0,0) with 
p = 0.6 in the P-AP configuration. 

Figure 3(b) exhibits the results for the transport regime (1,0). It is obvious that the 
spin-resolved CCs are all positive due to dynamical channel anti-blockade; moreover, the 
dynamical spin blockade strongly enhances S^^O) and weakly suppresses Stj.,H(0) and 
<Sti,6T (0) in comparison with the results of the paramagnetic case. Therefore, the two effects 
are positively additive and they result in an enhanced positive CC for the transport regimes, 
(1,0) and (0,1). For the double occupation case (1,1), inter-dot Coulomb blockade is the 
dominant mechanism governing transport (over the dynamical spin blockade), giving rise to 
a negative value, but enhanced in magnitude, for the up-up CC [Fig. 3(c)]. 

On the contrary, the situation is quite different when one of the circuits is in the AP 
configuration. For the AP case, the fast electron has spin up at the input terminal but 
has spin down at the output terminal, leading to an accumulation of up spins and a rather 
weak bunching of tunneling events associated with down spin at the output terminal. When 
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the top circuit is in the P configuration while the bottom circuit is in the AP configuration 
(P-AP configuration), we find very small positive values for 5V|-,ftj(0) and •S'tT,fai (0) , and thus 
a vanishing Sf&(0) for the transport regime (0,0) [inset of Fig. 3(a) and Fig. 2(d)]. A similar 
result is obtained when both circuits are in the AP configuration [Fig. 2(c)]. 

Our results show that a parallel-coupled QD in the P-P configuration is optimal for the 
occurrence of positive current CC: in the Coulomb anti-blockade regimes, (1,0) and (0,1), 
the maximum value of the CC can be nearly two times higher than that of the paramagnetic 
case. For instance, Sf&(0) ~ O.O67 with p = 0, while S t b(0) ~ O.I87 with p = 0.6 and ~ O.I7 
with p = 0.4. In the Coulomb blockade regime, (0,0), the current CC can be S t b(0) ~ O.O87 
for p = 0.6 and ~ O.O27 for p = 0.5. The sign-change of the current CC between different 
output terminals due to dynamical spin blockade has not yet been observed experimentally 
for the original setup, a three-terminal QD with ferromagnetic leads. 6 The main difficulty 
may be due to attaching the third ferromagnetic lead to the QD such that its tunnel coupling 
strength is comparable to that of the others. Our present proposal involves only a simpler 
element, an interacting QD connected to two ferromagnetic leads, which has already been 
realized experimentally— and it may provide a more accessible setup for the observation of 
positive CCs originating from dynamical spin blockade. 

In conclusion, we have analyzed the gate- voltage dependent CC noise between two output 
terminals in tunneling through two capacitively coupled QDs connected to four independent 
ferromagnetic electrodes. We find a sign reversal of the zero-frequency CC noise in the 
dynamical channel blockade regime, if the polarization of the electrodes is sufficiently strong. 
Moreover, in the dynamical channel anti-blockade regimes, positive CC noises are obviously 
enhanced in comparison with the results of paramagnetic leads. These results may be 
ascribed to the joint effect of dynamical Coulomb blockade and dynamical spin blockade, 
and are found to be quite robust against the spin-flip relaxation. 
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